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Abstract 

This paper presents closed-form analytic solutions to two illustrative problems in solar physics that have 

been considered not solvable in this way previously. Both the outflow speed and the mass loss rate of the 

^ . solar wind of plasma particles ejected by the Sun are derived analytically for certain illustrative approxi- 

o, 

CN , mations. The calculated radial dependence of the flow speed applies to both Parker's isothermal solar wind 

^ ' equation and Bondi's equation of spherical accretion. These problems involve the solution of transcendental 
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equations containing products of variables and their logarithms. Such equations appear in many fields of 
physics and are solvable by use of the Lambert W function, which is briefly described. This paper is an 
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\^ ' example of how new functions can be applied to existing problems. 
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I. INTRODUCTION 

Most stars eject matter from their atmospheres and fill the surrounding space with hot, low- 
density gas.'^ Astronomers and space physicists have studied the continuously expanding solar 
wind of charged particles from the Sun for almost a half century.^'^"^ The study of the solar wind 
as a unique plasma laboratory is compelling for several reasons. Many basic processes in various 
fields of physics (for example, plasma physics, electromagnetic wave theory, and nonequilibrium 
thermodynamics) have been detected in the solar wind and almost nowhere else. The solar wind is 
the closest example of a stellar wind, and stellar winds affect the long-term evolution of galaxies 
by injecting large amounts of matter and energy into the interstellar medium. On the more practical 
side, when solar wind particles impact the Earth's magnetosphere, they can interrupt communi- 
cations, threaten satellites and the safety of orbiting astronauts, and disrupt ground-based power 
grids. ^ 

The crown-like solar corona seen during a total eclipse is the place where the solar wind un- 
dergoes its initial acceleration. The shimmering auroras seen in northern and southern skies are 
the end products of the interaction between incoming solar wind particles and the Earth's mag- 
netic field. Sightings of the corona and the aurora go back into antiquity, but the first scientific 
understanding of the solar wind came at the beginning of the 20th century. Researchers gradually 
realized that there were strong correlations between the appearances of sunspot activity, geomag- 
netic storms, auroras, and motions in comet tails. In 1958, Eugene Parker^-^ combined these 
empirical clues with the knowledge that the bright solar corona consists of extremely hot (10^ K) 
plasma and postulated a model of a steady-state outward expansion from the Sun. Parker's key 
insight was that the high temperature of the coronal plasma provides enough energy per particle to 
overcome gravity and produce a natural transition from a subsonic (bound, negative total energy) 
state near the Sun to a supersonic (outflowing, positive total energy) state in interplanetary space. 
This theory was controversial at the time, but Parker had only to wait four years until the existence 
of the continuous, supersonic solar wind was verified by the Mariner 2 probe in 1962.^ 

Over the past decade, our understanding of the physics of the solar wind has increased dra- 
matically from both new space-based observations and the rapid growth of computer power for 
simulations (see Refs.^9,10 for recent reviews). The Ulysses spacecraft, for example, was the 
first probe to venture far from the ecliptic plane and soar over the solar poles to measure the solar 
wind in three dimensions r^ Remote observations of the solar corona have become significantly 
more detailed with data from space-based telescopes pouring in as never before. Figure 1 shows 
a snapshot of the corona as observed in 1996 by two instruments on the SOHO (Solar and Helio- 
spheric Observatory) spacecraft. However, progress in solar wind research often requires substan- 
tial numerical analysis, because even the most basic problems have not been tractable by analytic 
means. This paper takes advantage of a new transcendental function that is unfamiliar to many 
physicists — the Lambert W function — to illustrate how two fundamental solar wind problems can 
be solved analytically. 

In this era of efficient numerical computation, it is worthwhile to list the various ways that 
analytic solutions for the properties of the solar wind can be useful (beyond their pure aesthetic 
appeal). Analytic expressions often are used as initial guesses for more complicated iterative, time- 



dependent, or multi-dimensional calculations. Closed-form solutions also make it easier to study 
linearized perturbations to a known background state. The rapid evaluation of a large number of 
cases is facilitated by having analytic formulae, especially because many symbolic computation 
packages already contain optimized routines for the Lambert W function. Finally, the ability to 
write down simple expressions for solar wind plasma properties may make the extrapolation to 
other stars more tractable and physically understandable. 

This paper presents a brief overview of the Lambert W function in Sec. II and a summary 
of the governing equations of the solar wind in Sec. III. The use of the W function in solving 
the classical Parker solar wind problem, that is, the radial dependence of the wind speed for an 
isothermal plasma, is given in Sec. IV. The use of this function in solving for the mass loss rate of 
the solar wind is given in Sec. V. Conclusions and other potential applications of this function are 
given in Sec. VI. 



II. SOME PROPERTIES OF THE LAMBERT W FUNCTION 

Like many mathematical functions, the Lambert W function was derived and used indepen- 
dently by several researchers before the mathematics and computer science community settled on 
a common notation in the mid- 1 990s. '^ This function has been used to solve problems in elec- 
trostatics, statistical mechanics, general relativity, radiative transfer, quantum chromodynamics, 
combinatorial number theory, fuel consumption, and population growth ,'^i^?i^'^i'^i^^ but is still not 
widely known by physicists. 

The Lambert W function is defined as the multivalued inverse of the function xe^. Equivalently, 
the multiple branches of W are the multiple roots of the equation 

W(z)e^'^''^ = z, (1) 

where z is in general complex. There are an infinite number of solution branches, labeled by 
convention by an integer subscript: Wk(z), for fc = 0, ±1, ±2, ... If z is a real number x, the only 
two branches that take on real values are Wo(x) and W-i(x). Figure 2 plots these two branches, 
which are the only ones that are needed in the applications of this paper. 

Numerous formulae for the differentiation, integration, and series expansion of W are given in 
the references cited above (for example, Ref. Il2lll4r) . One useful result, which is applied in the 



following, is given here. Near the branch cut point at x = — 1/e, Wq = W-i = — 1, and the two 
real branches can be approximated to lowest order by 

Woix) ^ -X + ^l^lex (2) 

W-xix) ^ -1 - V2 -1- lex. (3) 

A useful way for expressing the solutions to a standard family of transcendental equations in 
terms of the Lambert W function is to note that the equation '^ 

\a{A + Bx) + Cx = \aD, (4) 



where A, B, C, and D do not depend on x, has the exact solution 
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The choice of solution branch usually depends on physical arguments or boundary conditions. 

III. GOVERNING EQUATIONS OF THE SOLAR WIND 

The expansion of the solar wind is treated traditionally as a problem of steady-state hydro- 
dynamics. The Sun has a strong magnetic field that confines and directs the flow of plasma (see 
Fig. 1), but along the magnetic "flux tubes," the dynamics is essentially independent of the strength 
of the field. Thus, the overall properties of the solar wind can be determined by solving the hydro- 
dynamic equations of mass, momentum, and energy conservation. ^-^^^ 

The equation of mass conservation for a single-component fluid is given by 

|^-hV-(pu) = 0, (6) 

where p is the mass density and u is the velocity. The equation of momentum conservation is 

^-h(u-V)u = — VP-hg, (7) 

at p 

where P is the gas pressure and g is the net external force on a parcel of gas, here assumed to be 
only due to gravity. The equation of total energy conservation, 

contains several terms that arise from the ionized nature of the near-Sun plasma. The solar corona 
exhibits a temperature of about 10^ K, which is at least two orders of magnitude higher than the 
temperature at the base of the atmosphere (that is, the photosphere and chromosphere). A major 
unsolved problem of solar physics is to explain what physical processes lead to such large amounts 
of energy into the coronal plasma. Significant progress has been made, though, by constraining 
the input energy flux density Fh empirically, even though the exact physical origin of this heating 
is not known. 

Some of the energy deposited into the corona is transported downward to the lower atmosphere 
via heat conduction (that is, the radial component of F^ is negative), and some of it is converted 
back and forth between kinetic energy, thermal energy, and gravitational potential energy. (See 
the final three terms in square brackets in Eq. dH), where G is the gravitational constant and Mq 
is the mass of the Sun.) Some of this energy also is lost in the form of radiation, because many 
of the free electrons that exist in an ionized plasma undergo collisions with bound atoms and lib- 
erate a fraction of their energy to photon emission. The radiative loss function $(T) encapsulates 
the elemental composition and atomic physics of the radiating coronal plasma as a function of 
temperature T. 



Equations ©-(El) can be simplified in several ways without sacrificing realism. We can neglect 
the partial time derivatives and restrict our analysis to time-steady solutions of a continuously 
expanding solar wind. Also, the vector terms can be expressed in spherical polar coordinates, as- 
suming that the variations exist only along the radial direction r and that all vectors have nonzero 
components only in this direction. For example, the radial component of the gravitational acceler- 
ation g is simply g,. = —GMq/t'^. 

Very near the Sun, the assumption of radial flow is not accurate because of the complex multi- 
pole structure of the solar magnetic field (see Fig. 1). At larger distances, though, the radial flow of 
the solar wind has been largely confirmed by spacecraft measurements.^-^ However, the Sun's slow 
rotation (once every 27 days) causes the flow direction and the magnetic field direction to become 
misaligned in interplanetary space. Because of the high conductivity of the solar wind plasma, the 
magnetic field lines become "frozen in" to the flow, that is, the magnetic field becomes a passive 
tracer of the flow, like drops of ink in a flow of water. A possible outdoor demonstration of this 
effect can be performed with a persistent, rotating source of water flow, like a lawn sprinkler. De- 
spite the fact that all of the water droplets are flowing radially away from the center, a snapshot 
at any time shows them arranged in a spiral "streakline." This field is completely analogous to 
the Parker spiral magnetic field pattern in the solar wind, which carries the imprint of the Sun's 
rotation, but still channels the particle flow to be radial. 

For the useful assumption of radial flow, the mass conservation equation is expressed in spher- 
ical symmetry as 

i-^(pMr2) = 0. (9) 

Because the radial derivative of pur^ is zero, this quantity is constant. We thus define the total 
mass loss rate from the entire Sun (in units of kg s^^) as M = Anpur^, where the factor of 4% 
comes from integrating over the full solid angle of the spherical Sun. The energy conservation 
equation is written as 

i d ( .. ^ ,u^ 5P GMq,.^ . _ 

^^,{AF«^Fc.pa(-.---Sl)]} = -/HT). (10) 

with the radial components of vectors written as scalars with the same notation. 

The gas pressure P can be eliminated from the momentum equation by applying the ideal gas 
law, assuming the fluid consists of a single particle species with mass m, 

P = ^^pa\ (11) 

m 

where k is Boltzmann's constant and a is an effective sound speed. For a hydrogen plasma, m is 
essentially the proton mass. An added assumption that simplifies the subsequent analysis is that the 
hot corona is isothermal, that is, that after the coronal heating takes hold, the ~10^ K temperature 
remains roughly constant as a function of radius. It is known from spacecraft measurements that 
the plasma temperature drops only by a factor of 10 from the inner corona (r ^ 1.5 -R0) to the orbit 
of the Earth (r ^ 215 Rq = 1 AU), where Rq is the solar radius.- Therefore, for the acceleration 
region of the solar wind (1.5 to IORq), the isothermal approximation seems sufficiently valid. 



If we substitute these conditions into Eq. Q and use the mass conservation equation, we obtain 



a^.du _2a^ GMq 
u dr r r^ 



, a .au La KjIVIq ^^_ 



It is noteworthy that the momentum conservation equation is now a true equation of motion be- 
cause the mass density p no longer appears. 



IV. THE PARKER SOLAR WIND PROBLEM 

The fluid in a steady-state stellar wind accelerates from rest to an asymptotic "coasting" speed 
far from the star, where the star's gravity has become negligible. This situation can only be main- 
tained by a gradual transition from a hydrostatic force balance close to the star (that is, where 
inward and outward forces cancel) to a net outward force at larger distances.'^ Parker^ recognized 
that this transition occurs naturally for a hot (million K) corona, where the gradient of the large 
gas pressure plays the role of the increasing outward force. Equation (IT2b shows the primary man- 
ifestation of the gas pressure gradient as the first term on the right-hand side, which for a constant 
a eventually must overtake the more steeply decreasing gravity term and result in a net positive 
(outward) acceleration. Interestingly, the dynamics described by Eq. (IT^ does not depend on how 
the corona is heated, but merely on the fact that it is heated. 

Parker also noticed that Eq. (IT2b exhibits a potential singularity at the "sonic point," u = a, 
because when this condition applies, the term in parentheses on the left-hand side is zero, and for 
an arbitrary radius (that is, a finite value for the right-hand side) the first derivative of the velocity 
du/dr must be infinite. However, there is one specific value for r where the right-hand side is zero 
as well. If the sonic point occurs at the critical radius Vc = GMQ/(2a^), then du/dr may remain 
finite and the wind solution remains physically realistic. Mathematically, this solution represents 
an X-type singular point, at which two solution trajectories in (r, u) space intersect with slopes 
of opposite signs and other solutions are hyperbolic about this point. The joint set of conditions 
r = Tc and u = a often is called the Parker critical point. 

There are two possible solutions that pass through the critical point: one representing a continu- 
ously accelerating outward flow of gas (the wind), and one representing an outwardly decelerating, 
but inward flow of gas (steady spherical accretion). Parker's wind solution was criticized initially 
for being too "finely tuned" because it seemed unlikely that a wind would naturally want to ac- 
celerate through the sonic point exactly at r = r^. However, it has been noticed recently that 
Parker's critical solution is the only truly stable wind solution to Eq. (IT2t . and all other outwardly 
flowing solutions are unstable. ^° Note that six years before Parker, Bondi^' recognized that the 
inward accretion solution also could represent real astrophysical flows. It was also found that this 
solution, like Parker's, is a stable attractor and represents the maximum amount of mass that can 
be consumed (in steady state) from an external source.— 

Equation (IT2ll is a first-order ordinary differential equation that is separable. The integration of 
the left side from a to an arbitrary u and the right side from r^, to an arbitrary r yields an implicit 



transcendental equation for u and r: 

?y T / 1 1 \ 
{u^ - a^) - a^ In (— ) = Aa^ In (-) + 2GMq ( ) . (13) 

We rearrange terms and define tlie dimensionless variable y = {u/af, so that Eq. (IT3t becomes 

\ny -y = \nD{r), (14) 

where 

I^(r)=(-)"'exp[4(l-^)-l]. (15) 

Tc T 

Thus, using Eqs. (JH) and (O, the Parker/Bondi solutions have the general analytic solution y = 
—W[—D{r)]. For all values of r, D{r) ranges between and 1/e, so we must choose between the 
two branches Wq and W^i (see Fig. 2). Additionally, the full solution involves one choice below 
the critical point and the opposite choice above it. The accelerating Parker solar wind solution is 
given specifically by 

, \ -a^Wo[-Dir)] r <r, 
u^ = I (16) 

\-a^W-i[-Dir)] r>rc 

and the opposite choices must be made to obtain the Bondi accretion solution. 

Figure 3 shows a set of solutions for the Parker solar wind with six choices for the constant 
sound speed a. These solutions are compared to curves showing empirical (that is, observationally 
derived) speeds for the fastest and slowest types of solar wind flow that have been seen. Because 
our only direct measurements of the wind speed have been exterior to the orbit of Mercury (r > 
60Rq), indirect methods are needed to determine the wind speed at distances closer to the Sun. The 
analysis of ultraviolet photons emitted from the corona has provided new ways of probing the solar 
wind's acceleration, '° but a more traditional method is to measure the density of particles in the 
corona and use mass conservation, that is, the steady state version of Eq. Q, to compute the wind 
speed. The density can be measured by observing the linear polarization of Thomson-scattered 
visible light in the corona; the degree of polarization is directly proportional to the number of free 
electrons along the line of sight. 

In the following we give a simple parameterization^^ of the radial dependence of density as 
observed in the source regions of the fast and slow components of the solar wind (at the minimum 
of the Sun's 1 1-year magnetic cycle): 

p (fast) ^ 2.37 X 10-1^ (1 + ^ + ^) g/cm^ (17a) 

p(slow) ^ 6.21 X 10-1^^ + ^ + ^) g/cm^, (17b) 

\ rfiZ, rfj /-y»D ' 

iXj iXj JU 

where x = r / Rq. The wind speed at any radius is thus proportional to p~^r~^ times a normal- 
ization constant that is given by specifying the measured wind speed at 1 AU. Note that at large 
distances the dominant terms in Eq. (llTt will be the l/x^ terms, and thus u oc p^W~^ at large 
distances will approach a constant coasting speed. 
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The density is generally higher in the slower component of the solar wind, which emerges 
mainly from bright "streamers" around the solar equator and reaches speeds at 1 AU of 300 to 
500km/s. The density is lowest in the fast solar wind that emerges mainly from dark "coronal 
holes" at the north and south poles and reaches speeds at 1 AU of 600 to 800km/s.'' Figure 3 
shows that the acceleration of the slow wind has a very similar shape to the analytic solutions 
given in Eq. STE\i . The fast wind has a slightly steeper profile in the corona because this plasma 
is not isothermal, and because it also flows slightly "super-radially" (that is, the magnetic field 
over the poles flares out like a trumpet and the equations are not represented exactly by spherical 
symmetry; see Fig. 1). However, much of the essential physics of solar wind acceleration remains 
encapsulated in the radial, isothermal problem. 

One practical benefit of having an analytic expression for u{r) is being able to easily find 

asymptotic expansions for various limiting cases. In the nearby vicinity of the critical point, that 

is, for 1 1 — (rc/r)| <C 1, the series expansions given by Eqs. Q and © can be used in conjunction 

with the series expansion of D{r) about the critical point to obtain a single expression for radii 

near Vc'. 

I 2r^ 
M~ay3 -, when r ^ Tc. (18) 

Other expansions can be used to obtain approximations for r <^ r^ and r ^ Vc- 



V. THE MASS FLUX PROBLEM 

The above solution (Eq. lfT6l ) for the solar wind speed u{r) is only half of the problem. Because 
the mass density p was eliminated from the equation of motion, we know how fast the gas is 
accelerating, but we do not know how much gas is being ejected. The determination of the solar 
wind mass loss rate M is the second half of the problem which, interestingly, also is addressable 
using the Lambert W function. 

The Sun is observed to lose mass at a rate of approximately 10^'"* solar masses per year 
(M0/yr). This unconventional unit is useful because it can be compared easily to a firm upper 
limit derivable by dividing the mass of the star by its lifetime. For the Sun, with an expected 
main-sequence lifetime of about 10^° years, this upper limit is of order lO~^°M0/yr. Thus, the 
solar wind is expected to drain away no more than one ten-thousandth of the Sun's mass over the 
next few billion years. (Some hotter stars lose mass at much higher rates, with the wind having a 
substantial impact on the star's late stages of evolution.-) 

There is still not universal agreement about what determines the Sun's mass loss rate. ^'^i^^i^^i^^ 
Our analytic solutions apply to only one of the several suggested mechanisms. In this class of 
radiative energy balance models, first outlined in detail by Hammer,— M is determined at the base 
of the corona by the interplay between the heating and cooling terms in Eq. (ITOb . the equation of 
energy conservation. Because we have solved for u(r), the determination of the mass loss rate 
requires only the solution for the density at a single radius. 

To simplify Eq. (ITOb further, the solar atmosphere can be considered to consist of two concentric 



layers: the cool (^10"^ K), high density chromosphere, and the overlying hot (~10^ K), low density 
corona. The transition between these layers has been observed to be exceedingly thin — about 0.1% 
of a solar radius — so that the radial derivative in Eq. (ITOI) can be expressed as a simple difference 
of quantities above and below the transition zone. Because of the relative thinness of this zone, 
we can ignore both the small change in the gravitational potential energy between the two layers 
and the spherical divergence, that is, the r^ terms inside and outside the braces. Also, the kinetic 
energy term u^ /2 can be ignored because the solar wind speed has been seen to be negligibly 
small (that is, very subsonic) at the solar surface. Finally, the coronal heating term itself, Fh, is 
ignorable because we are concerned with layers below where the majority of the heat is deposited. 
Thus there are only three dominant terms in the energy balance: 

^{Fc + SnukT) = -p^^{T), (19) 

dr 

where for convenience we rewrite the mass density p as the product of the particle mass m and a 
number density n, that is, number of particles per unit volume. In summary, at the coronal base 
the heat is conducted downward from where it is initially deposited, some of it resides at the base 
as enthalpy, and the remainder is lost as radiation. 

The steady state balance of mass and momentum across the thin transition zone also de- 
mands that the products nu (mass flux) and nT (gas pressure) remain roughly constant. This 
mass flux constraint is used, together with an empirical form^^-^^- for the radiative loss function 
p2$ = n^AT-^/^ (where A=1.9x 10"^^ Wm^ K^/\ to obtain 

^ + 5nuk^ = -An^T-'/\ (20) 

dr dr 

The differential equation d^Ut is transformed by multiplying both sides by the heat conductive flux 
Fc- Note, though, that it is advantageous to multiply the left-hand side by Fc itself and to multiply 
the right-hand side by the definition of the classical conductive flux, 

Fc ^ -«:TV2^, (21) 

dr 

where k is the Spitzer-Harm^° heat conductivity in an ionized plasma, which has a value of 8.8 x 
10^'^ Wm^^ K^^/^ for the range of densities and temperatures of the corona. We rearrange and 
divide all terms by a factor of dT/dr and obtain the following form of the energy balance equation: 

dF 
iFc + Fc^ = ^, (22) 

where the quantities ^ = 5nuk and ^ = n^nAT^ are assumed to be constant across the thin 
transition zone. 

The above form (Eq. |I23) of the energy equation is separable and integrable with T and Fc as 
the independent and dependent variables, respectively. Once integrated across the transition zone, 
though, the full equation contains terms for T and Fc in both the upper and lower layers. The 
terms corresponding to the lower (chromospheric) layer can be neglected because the values of 
both T and Fc are several orders of magnitude smaller in comparison to their counterparts in the 



upper (coronal) layer. Thus, the integrated transcendental equation relates the values of T and Fe- 
at the coronal base to one another and is independent of their values in the chromosphere: 

Fc + ^T+^\n{l-^)=0. (23) 

We note that both ^ and ^ contain the number density at the coronal base (which we call no). 
Hence, the Lambert W function can be used to solve for this quantity, with 

5uokFc 

and the argument u of the W function is 

,= -exp(-?^)-l. (25) 

In this formulation, the only "free" variables are Uq, Fc, and T (all evaluated at the coronal 
base). Just as in the Parker solar wind application, the argument cu falls between — 1/e and 0, thus 
making the choice between the Wq and W^i branches necessary. In this case, though, the physical 
choice (VT-i) is apparent because the Wo branch gives a negative density. 



To use the above solution (Eq. lEH) to compute the mass loss rate, we use the results of SecHVl 
to fix Mo for a given isothermal Parker wind model. A realistic median value of T is 1.2 x 10^ K, 
which has an outflow speed uq of 0.96km/s at r = Rq. To estimate the applicable values of Fc at 
the coronal base, Eq. (OTT) can be solved assuming a finite-difference temperature gradient across 
the thin transition zone. A thickness Ar of 0.001 Rq gives values of Fc of order -1000 Wm~^. 
In more accurate models,^*^ though, the temperature gradient is a bit less steep at the top of the 
transition zone, and Fc ranges between -50 and -200 Wm^^. For this range, Eq. (l24b for no 
yields values of the number density between 6 x 10'^ and 3 x lO^'* m^^. The mass flux, integrated 
around the whole sphere, is then M = ATrmpTioUoRQ, and the resulting values range between 
9 X 10^^^ Mq/jv and 4 x 10^^'^ Mq/jt. The observed solar mass loss rate is observed to vary 
between about 2 x 10^"^ M^lyx at the solar minimum and a few times that at the solar maximum. 
If we take into consideration the large number of approximations we have applied, we conclude 
that the agreement is good. 



VI. CONCLUSIONS 

We have presented new analytic solutions to two simple problems in solar wind physics. The 
Lambert W function used in these solutions was defined and publicized only about a decade ago, 
but it has rapidly become a convenient tool for mathematical physicists. The elegance of explicit 
solutions to equations thought previously to be expressible only implicitly is clear, but there also 
are many practical benefits to having explicit solutions as well. 

There are other potential applications of the Lambert W function in solar and space physics. 
A transcendental equation solvable in terms of W arises in a calculation of the electric potential 

10 



drop that exists between the Sun and the edge of the solar system.^^ (An ionized plasma exhibits 
local charge neutrality because of electrostatic screening, but for solar wind particles in the Sun's 
gravitational field this neutrality is possible only by setting up a radially varying electric field.) 
Functions with temperatures T appearing both inside and outside exponents occur when calcu- 
lating the energy distributions of solar photons (for example, the Planck blackbody function) and 
electrons in excited atoms (for example, the Saha ionization equation). The Lambert W function 
can thus be used in a variety of ways when the need to solve for T arises. Further applications are 
expected to clarify the physics of many types of systems. 
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Figure Captions 




FIG. 1: The solar corona on August 17, 1996, with bright regions plotted as dark. The inner image was taken 
by the BIT (Extreme-ultraviolet Imaging Telescope) instrument on Solar and Heliospheric Observatory 
(SOHO), and is sensitive to the ultraviolet emission of Fe^^' ions at temperatures of about 10^ K. The outer 
image was taken by the UVCS (Ultraviolet Coronagraph Spectrometer) instrument on SOHO by blocking 
out the bright disk to see the much dimmer ultraviolet emission of O"*"^ ions at temperatures exceeding 10^ K. 
The magnetic field lines are from a model of the corona at the minimum of its 1 1-year activity cycle.— 
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FIG. 2: The two real branches of the Lambert W function. 
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FIG. 3: Analytic solutions of the isothermal Parker solar wind equation, plotted as outflow speeds versus 
heliocentric distance (in units of a solar radius; Rq = 6.96 x 10^ m). Individual solutions (thin solid lines) 
are labeled with the locations of the critical point (circles). From bottom to top, the modeled coronal tem- 
peratures are 0.1, 0.2, 0.5, 1, 2, 5MK, respectively. Shown for comparison are observationally constrained 
wind speeds for polar (thick dashed line) and equatorial (thick dotted line) flow at solar minimum. 
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